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Abstract 

On the basis of the closed-time path formalism of non-equilibrium quantum field theory, we derive 

the real-time quantum dynamics of heavy quark systems. Even though our primary goal is the 

description of heavy quarkonia, our method allows a unified description of the propagation of single 

| heavy quarks as well as their bound states. To make calculations tractable, we deploy leading- 

P-l 

order perturbation theory and consider the non-relativistic limit. Various dynamical equations, 

a? 

J3 such as the master equation for quantum Brownian motion and time-evolution equation for heavy 

quark and quarkonium forward correlators, are obtained from a single operator, the renormalized 
effective Hamiltonian. We are thus able to reproduce previous results of perturbative calculations 
of the drag force and the complex potential simultaneously. In addition, we present stochastic 
time-evolution equations for heavy quark and quarkonium wave function, which are equivalent to 
the dynamical equations. 

PACS numbers: 
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I. INTRODUCTION 



Vigorous study of the quark-gluon plasma (QGP) has revealed its strongly coupled nature 
This novel state of matter, composed of deconfined light quarks (up, down, and strange 
quarks) and gluons, which once existed around 10 -4 to 10 -5 second after the Big Bang, is 
expected to be recreated at the Relativistic Heavy Ion Collider (RHIC) and at the Large 
Hadron Collider (LHC). Relativistic hydrodynamics has been used to describe the space- 
time dynamics of the bulk matter in such collisions j^j, which when compared to elliptic 
flow data from RHIC revealed a very small shear viscosity (rj) to entropy density (s) ratio 
rj/s ~ (1 - 2.5)/47r {3]. The dimensionless ratio r\js has since become one of the most 
sensitive measures of strongly coupled nature of QGP. In addition to the elliptic flow of 
light hadrons, which is sensitive to the collective behavior and the initial geometry, other 
observables such as electromagnetic radiation, jets, and heavy quarks also provide valuable 
insights to the property of QGP. Among these probes, we concentrate on the heavy quark 
probe of QGP. 

In relativistic heavy-ion collisions, heavy quarks (charm and bottom quarks) are pro- 
duced in initial hard partonic processes and are expected to probe the properties of QGP 
through their subsequent evolution in the medium. Better theoretical understanding of how 
the interaction of heavy quarks and medium particles modifies the kinematics and dynam- 
ics of the heavy quarks inside the strongly coupled QGP fireball will allow us to utilize 
them as well controlled hard probes (4]. Heavy quarks have been considered to probe trans- 
port and static screening properties of the medium. Momentum spectra of heavy quarks 
are modified due to energy loss and diffusion processes in the medium Survival of 

heavy quark-antiquark bound states (J/^f and T), or heavy quarkonia, is expected to be 
sensitive to how the potential is altered by the Debye screening of the color charges [7, 8], 
although the recombination of heavy quarks into bound states might compensate the pos- 
sible suppression of heavy quarkonia |9|, llO[. Experimental data suggest that all of these 
phenomena actually take place in heavy-ion collisions at LHC and RHIC. The suppression 
of the nuclear modification factor and elliptic flow of open heavy quark observables, namely 
single electrons, single muons, and D mesons, indicates a strong drag force acting on the 



heavy quarks 



33 



Sequential suppression of T(1S_ L 2S_ L 3S) is observed |37H39j|. while the 



recombination processes may be important for J/ty 40H44j . This suggests the usefulness of 
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the T family as a clean thermometer of the QGP and also the necessity of knowing wave 
functions of charm quarks at hadronization. Motivated by these phenomenological inter- 
ests, theoretical computations have been performed and are under active development for 



the diffusion constant Jl 
for current correlator 125 



-|22] , heavy quark-antiquark free energy 23_|, |24[ , spectral 



27] . and complex potenti al I28H32I ]. Of these, some 



are based on perturbation theory, some 
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19 



23 



27 
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unction 
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others 



20-22 



31] on lattice QCD simulation, and 



32] on the conjecture underlying the AdS/CFT correspondence. 



In this paper, we will present a new derivation of the dynamical description of heavy 
quarkonia in a thermal medium, which might also lead to a better description of the re- 
combination process of heavy quarks. To discuss the dynamics of heavy quarkonia in the 
medium and recombination of heavy quarks into bound states, a quantum description is 
indispensable. The most convenient framework to study such problems, namely the real- 
time dynamics of a q uantum system interacting with a thermal environment, is that of open 



quantum system 



45 1 . Presently, there exist a few publications based on this approach in 



similar contexts 46-50]. In Ref. [46, 147(| , analyses using a QCD version of the quantum opti- 
cal master equation are presented. It is important to note that the quantum optical master 
equation is derived under the rotating wave approximation, which assumes the following 
separation of time scales te, t$ <C tr 45] . Here te denotes the correlation time for the 
environment, rs ~ — u'\ with u and u' being the typical frequencies or energy levels of 
the open quantum system, and tr the relaxation time for the open quantum system. Since 
the relaxation time of the system scales as tr oc M/T 2 , this approach should be ideal in 
studying the real-time dynamics involving transitions among well-separated bound states. 
On the other hand, quantum Brownian motion, which is studied in Refs. 48|, |49J, is derived 
by assuming a different time scale separation te <C ts, tr 45j. The spacing of energy levels 



of the system is expected to get smaller when the heavy quarks are almost liberated from 
bound states, be it in a heavy-light meson or heavy quarkonium, where thus this approach 
will become suitable. Note that quantum Brownian motion of heavy quarks should not 
necessarily be described by the Caldeira-Leggett model 51|. 



We will follow and extend the philosophy of the stochastic potential 



50] , which has several 



appealing features. The imaginary part of the complex heavy quark-antiquark potential 
3l] can be naturally interpreted as damping due to wave function decoherence between 



different realizations of the stochastic potential. It models the random microscopic collisions 
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between heavy quarks and medium particles and thus makes the energy of the heavy quark 
system fluctuate. In each realization of the stochastic potential, the particle number of heavy 
quarks and that of heavy antiquarks are explicitly conserved separately, as they should 



be in the non-relativistic limit. In Ref. 



50], however, it is pointed out that the energy 



increases linearly without bound in the stochastic description and thus the system does 
not get thermalized. What is missing in the stochastic description 50( is an irreversible 
process, namely friction, which prevents the energy to rise forever. Therefore when the 
heavy quark mass is finite M < oo, the proper quantum description of heavy quark systems 
must incorporate this friction process as well as the thermal fluctuations, which are already 
properly modeled by the stochastic potential in the M — > oo limit. Interestingly, even 
though a classical description is mostly sufficient and consistent in describing energy loss 
and diffusion processes of single heavy quarks, its quantum description becomes a requisite in 
describing the quantum dynamics of heavy quarkonia. A similar argument but in a different 



context is found in Refs. 



52 



We will study how friction should be incorporated in the quantum description of heavy 
quark systems from first principle on the basis of the closed-time path (CTP) formalism 



of non-equilibrium quantum field theory 



55] . We apply the CTP formalism to QCD and 



derive the real-time dynamics of heavy quarks as open quantum systems, in leading-order 
perturbation theory and in the non-relativistic limit. The essence is encapsulated in a 
renormalized effective Hamiltonian for the doubled degrees of freedom on each time path, 
which provides all the time-evolution equations of interest: the master equation of the 
reduced density matrix and the time-evolution equation of the forward correlator for an 
arbitrary number of heavy quarks and heavy antiquarks. It is found that thermal fluctuation 
(friction) derives from the two-point functions of gluon fields in leading order (next-to-leading 
order) in the low frequency limit and that the order in the low-frequency expansion of the 
two-point functions directly translates into the order in the non-relativistic expansion of 



ieayy_ 
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quark interaction terms. By ana 



14] and the complex potential 28 



yzing these time-evolution equations, the drag force 



29] in leading-order perturbation are reproduced 



simultaneously. Concentrating on around the diagonal parts of the reduced density matrix 
in the position representation, the master equation is reduced to that of the Caldeira-Leggett 



model 



51] . We will then show how friction can be incorporated in the stochastic time 



evolution. We will see that complex noise as well as real noise are necessary in the stochastic 
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description. The stochastic evolution of the wave function and that of its conjugate, or in 
other words the forward and backward propagation of the wave function, turn out to be 
independent, which is convincing because friction processes are irreversible. 

This paper is organized as follows: In Sec. [Ill we wm introduce the close-time path 
formalism of non-equilibrium quantum field theory and derive both the influence functional 
and the effective action in leading-order perturbation theory and in the non-relativistic limit. 
In Sec. PH we will derive an effective Hamiltonian corresponding to the effective action. The 
normal-ordered renormalized effective Hamiltonian is also calculated. In Sec. IIVI various 
real-time dynamical equations, namely the master equation for the reduced density matrix, 
the time-evolution equation of the forward correlator, and a stochastic Schrodinger equation, 
are derived. In Sec. |V] we conclude our study and give future prospects. Throughout this 
paper, we adopt natural units h = c = k-Q = 1 and all operators in Hilbert and Fock spaces 
are denoted by bold fonts. 



II. INFLUENCE FUNCTIONAL OF HEAVY QUARKS 

In this section, we will derive the influence action for heavy quarks with mass M in a 
medium with temperature T in the leading order of g and up to the order y/T/M in the 
heavy quark mass. 



A. Path integral on a closed-time path 



In studying real-time dynamics of quantum systems, one often encounters a situation 
where in the Heisenberg picture iV-point functions (N = n + m) of the form 



Tr 



T{(f) n (t n ) ■ ■ ■ 4>2{h)4>x{h)} PtotT '{Vm(«m) ' " ¥>2 {s 2 )lfl (si) } 



Tr 



T {(pm(s m ) ■ ■ ■ (^(^Vl^l)} T i<i>n{tn) ' ' ' </>2 (h) 01 (*l) } Pt 



:d 



are of interest. Here T(T) denotes the operation of taking the time-ordered (anti time- 
ordered) product and p t ot is the density matrix of the total system. In the path integral 
formalism, such correlators are obtained by defining the fields to be used in the functional 
integration on a closed-time path contour C labeled bys (—00 < s < 00) as illustrated 
in Fig. [1] 56H58||. Correspondingly, source fields in the generating functional must also be 



5 



Im[t] 



C 



3* 



Re[t] 



FIG. 1: Closed-time path C in complex time plane, parameterized by s. 
defined on the same contour: 



Z[r],r]] 



V[ip,ip,q,q,A](ip,q,A (s = -oo)\p to t\^,q,A (s = oo)) 



x exp 



% j ds I d 3 x{C QC D(ip^,q,q,A)+r]'ip + 'ipri} 



(2) 



where (77,7/) are sources for heavy quark fields, (V>,V0 are the heavy quark fields, (g, q) are 
light quark fields, A = A a ^ are gluon fields, and £qcd is the QCD Lagrangian density. 
Here functional integration of the ghost fields and the Faddeev-Popov term are omitted in 
order to simplify the notation. For later purposes, we only introduce source fields for heavy 
quarks. Using the single time representation, this path integral can be rewritten as 

Zfai, 771,772,772] = 



D[V>i,^i,gi,gi,Ai] J T>[ip 2 , ^2,^2, 92,^2] 
x(ip 1 ,q 1 ,A 1 (t = -oo)|p tot |^2,?2,^2 (t = -oo)) 
x exp 

x exp 



% j d 4 x{C QC D(ipi,ipi,qi,qi,A 1 )+r)i'ip 1 +'ipiri 1 } 
-i / d A x {£qcd(^2,V'2, 92,^2, A 2 ) + 772^2 + ^2772} 



(3) 

with boundary conditions at t = 00: fe(oo)7 ] = 7/^(00), [51(00)7°]^ = g 2 (oo), and 
Ai(oo) = A 2 (oo). 

In the following, we consider a density matrix in the Heisenberg picture (initial condition 
of a density matrix in the Schrodinger picture) of the form p to t = Pe ® PS) where p^? is 
the equilibrium density matrix for the environment degrees of freedom, i.e. light quarks and 
gluons, and pg is an arbitrary density matrix for the system degrees of freedom, i.e. heavy 
quarks. With this density matrix, the complex weight for the initial fields at t = —00 is 
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given by 



(V>l,<?l,Ai (-Oo)|p tot |V>2,<?2,^2 (-OO)) 

= (qi,Ai (-oo)\p c £\q 2 , A 2 (-00)) • (^i(-oo)|p s |^ 2 (-oo)). 



(4) 



B. Integrating out the environment degrees of freedom 

Since we are interested in the dynamics of heavy quarks as an open quantum system, let 
us integrate out the gluon and light quark fields: 

7?i, 773,7/2] = J P[^i,^i,^2,^2](^l(-00)|ps|^2(-00)) (5) 

xexp ij d 4 x {£^(^1,^1) + fjiipi + ^17/1} 

x exp -i J d 4 x {£^2, ih) + fj 2 i>2 + */> 2 %}J ■ 2fv$i*Wi, 

ZfvU^JT] = J ^[quqi,Ai,q 2 ,q2,A 2 \{q 1 ,A 1 (-oo)\ p c £\q 2 , A 2 (-00)) (6) 

x exp i J d 4 x {C g+q (q 1 , q^AJ- gj^A^} 

x exp -i J d 4 x {C g+g (q 2 , q 2 , A 2 ) - gf 2 A a 2p ) 

Here and C g+q are the Lagrangian densities of free heavy quarks and interacting glu- 
ons and light quarks. The couplings between heavy quarks and gluons appear in the 
Feynman- Vernon influence functional Z-pyiji^, j^] 59[ as source terms for gluons. Note 
that In ZfvL?!^, j 2 ^] is a generating functional for contour-ordered connected iV-point func- 
tions of gluons (N = n + m): 

1 \ n ( 1 \ 171 5 n 5 m 

ifj 



ig) Sjt 1 " 1 (x 1 ) ■ ■■5f 1 ^(x n ) 6j b 2 ^( yi ) ■ ■■5j b 2 ^{y m ) 
T [A b u \ ( Vl ) ■ ■ ■ A b Z (y m ) } T {A;\ ( Xl ) ■ ■ ■ A* (x n ) } 



1 ry 1 -an -au-i 

lnZ F vbi ,J 2 1 



(7) 



Then, if we assume (A aM (x)) = 0, the expansion of In Zpv [ji M 5 j/2 M ] in terms of and is 



lnZpvbT,^] 



9 - j d 4 xd 4 y(f/(x), j7(x)) 

Gl b ^{x-y) -G<^(x-y)\ ( j b »{y) 
-G^ix-y) Gi^(x-y) ) \j b 2 »(y) 

■ o r ,a M - Qui 



0(g 3 ) (8) 



where the (bare) two-point functions are defined by 

G F ab ,^(x -y) = (TAl{x)A h v {y)), G^x -y) = (T A^x)Al(y)), (9) 
G^(x -y) = (Al{x)A b M), G<^(x -y) = (A b u (y)A°(x)). (10) 

Although the argument so far is applicable to heavy fermions (muons) in QED straight- 
forwardly, there are a few differences to be remarked. In QED without light fermions (elec- 
trons), the expansion terminates at the second order in j± and j'2, while in QED with 
electrons, the expansion continues. Also, we need to subtract from jf _0 and j^ 1-0 the infi- 
nite background charge of heavy fermions which fill the Dirac sea. In QCD, the expansion 
continues even without light quarks, while heavy quarks in the Dirac sea do not contribute 
to the background color charge. 



C. Non-relativistic limit and long-time behavior 

The large separation of energy scales between the environment and the heavy quark 
system M ^> T allows us to treat heavy quarks as non-relativistic particles. In order to 



describe their dynamics, let us carry out the 



quark fields in momentum space ip{t,p) 



60 



6l|: 



dy-Wouthuysen transformation on the heavy 



with 7^ being the Dirac matrices in the Dirac representation. This transformation results 
in the following mass and kinetic terms expanded in terms of 1/M: 

= £™ + 0(T 2 /M), (12) 
Cf = Qt [ ido _ M + ^Q + Q c [id, + M-^j Ql (13) 

Here Q and Q c are the two-component Pauli spinors for a heavy quark and heavy antiquark, 
defined by ip' = T (Q, Ql). We make the parametric estimate VQ^j ~ VMT ■ Q^A since the 
typical momentum of heavy quarks is p ~ yMT in thermal equilibrium. In this approxi- 
mation, the kinetic term in CT 1 behaves as ~ (M + T). The non-relativistic limit of the 
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currents is in turn obtained as 



f° = Q ] t a Q + Q c t a Ql f = f v + f m (14) 

* = Q] eQ ~ Q < (%r) taQl + ° ( (T/M)3/2 ) ■ (15) 

f m = Q c at a Q + Qiat a Ql + O (T/M) . (16) 

Here a denotes the Pauli matrices. The vector meson term jm represents heavy quark 
pair creation and annihilation processes through the vector channel. We assume j a ~ j° 
which amounts to eliminating from the currents, equivalent to making the quenched 
approximation. The strength of the currents is then estimated as j a0 ~ 1 and j° ~ y/T/M. 
For later purposes, we also take the additional parametric estimate d j a0 = —V • j% ~ 



y^T/M, doj^ = 0, based on the free equation of motion 62 1. 

In order to study the heavy quark dynamics on long time scales, let us concentrate on the 
low frequency behavior of the gluonic two-point function (for a moment, the superscripts 
and subscripts are suppressed as far as possible): 

G(w,x-y) w G(x-y)+uG'(x-y), (17) 

/oo 
d(x°-y°)G(x-y), (18) 
-oo 



dtu 



-oo 

'■OO 



G\x-y) = ^-G(u,x-y) =i I d(x° - y°)(x° - y°)G(x - y). (19) 



w=0 



oo 



This is equivalent to the instantaneous approximation: 

G(x- y) « G(x - y)6(x° ~ V°) + iG'(x - y) J_ 5(x° - y°) , (20) 
with which the interaction term is approximated as 

d 4 xd 4 yj(x)G(x-y)j(y) (21) 
dt J d 3 xd 3 y G(x-y)j(t,x)j(t,y) - l -G'{x-y) {d j(t,x)j(t,y) - j(t,x)d j{t,y)} 
In covariant and Coulomb gauges, G^ = results from the tensor structures. Since we want 



to derive the influence action Spy [7°^, j^] at leading order in 0(g 2 ) and up to y/T/M, only 
Goo and G' 00 matter in these gauges. Note that we have made instantaneous approximation 
to describe the long-time behavior of heavy quarks and obtained the interaction term in the 
1/M expansion. 
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All of the gluonic two-point functions in the instantaneous approximation are expressed 
in terms of two real functions, G^ 00 (x — y) and G> b00 {x — y): 



G% {S-y) = G'! bfi0 (x-y) = 0, 
Gab,oo(% ~ v) — Ga b 00 (x — y), 
G'abw{x-y) ~- 



-iG% (x -y) + G> b n Q (x - y), 



x~y) = 2rpG> bQ0 (x-y), 



(22) 
(23) 
(24) 
(25) 



where we define the retarded propagator G^; b00 (x — y) = i9(x° — y°) ([Aq(x) , A^y)]). In the 



instantaneous approximation, G^ t 



afe.OO 



X 



y) and G> b00 (x — y) are expressed in terms of an 



Euclidean correlator and a spectral function: 



dr{A a0 (t = -ir,x)A b0 (t = 0,y, )), 



d 

r a b,oo(x-y) = T—a abfi o(uJ,x-y) 



oj=0 



dx ° e ^-y°) {[Aao{x ^ AbQm . 



Gab,00(x-y) 

G> 

cr abt00 (uj,x-y) 

The detailed derivation of these relations is given in the Appendix |A] 
In summary, we have derived an effective action S1+2 defined as 

5*1+2 = 'S'kin + SW 

= J d 4 x [£^(^1,^1) - £^(^2,^2)] + Spy, 
at leading order of the perturbative expansion and in the non-relativistic limit 

C CNR I CLONR 1 

^1+2 - ^ k i n + ^FV H 

d 4 x 



(26) 

(27) 
(28) 



CNR 
°kin 



CLONR 

°FV 



dtdVy [V(x - y)jf(t, x)jf(t, y) -V*(x- y)jf(t, x)jf(t, y)] 



(29) 

(30) 
(31) 



+ 1 I dtd 3 xd 3 yD(x-y)jf(t,x)jf(t,y) 



Mi 



■aOi 



where = Q 



1 



dtd 3 xd 3 yV x D(x - y) • |.# NR (t, x)jf(t, y) + jf(t, x) j 2 a ,NR^ y) , (32) 



aO/ 



■aOi 



2iM 



t a Q - Q c 



2iM 



t a Ql [67L Here we define a complex potential 



V(x — y) and a dissipative coupling D(x — y) through: 

-9 2 {G% (x-y) + iG> bt00 (x-y)} = V(x-y)5, 



ab i 



-9 G ab>00 (x - y) = D(x-y)S ab . 
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(33) 
(34) 



Note that the retarded propagator G^ boo (x — y) contributes both to the real part of the 
potential while the forward correlator G^ b00 (x — y) contributes to the imaginary part of the 
potential and the dissipative coupling. Since the heavy quark fields "0i,^i and ^2,^2 are 
coupled by D(x — y), we can interpret Eq. (I52"j) as a theory of two interacting species of 
heavy quarks. 



III. RENORMALIZED EFFECTIVE HAMILTON! AN 

In this section, we will construct a renormalized effective Hamiltonian from the effective 
action Eq. (130|) derived in the previous section. We use the notation ip = T (Q, Q\) for the 
4-component Dirac spinor after Foldy-Wouthuysen transformation. 



A. Effective Hamiltonian 



Here we derive the effective Hamiltonian H1+2 corresponding to the effective action 
5*1+2 = S^P + Spy NK + • • • for the interacting two species of heavy quarks. In this pro- 
cedure, we must take care of the time argument of each heav y q uark field since the opera- 
tors must be time-ordered back in the canonical formulation 63|. Although all the kinetic 



and interaction terms have the same time arguments, originally they were different. Time 
arguments in the bilinears forms of heavy quark fields, namely the kinetic terms and cur- 
rents should be read as ip\(t + e/2,x), ip\{t — e/2,x), ip\{t — e/2,x), ^(i + e/2,x) with 
e > 0. The products of currents should be understood as sums of all possible time order- 
ings with proper normalization, symbolically expressed as ^ S(Tes„{n" = ij(^cr(i), Xi)} with S n 
being the permutation group of n elements and a being one of these permutations, where 
all the time arguments are t ~ t\ > ti > ■ ■ ■ > t n ~ t. £py NR applies to n = 2 in this 
expression. The ordering in the fermion bilinear derives from the time-evolution of the den- 
sity matrix e~ ltHQCD p tot e ltHciCD in the Schrodinger picture and that in the current products 
derives from the current couplings which are nonlocal in time before making the instan- 
taneous approximation. In order to simplify the equations, we define a notation T sym by 

Now let us introduce new variables t/> 2 = or (<52,Q2c) = (Q^Qzc), which possess 
properties that will make proceeding arguments transparent. For example, the time argu- 



11 



ments in the fermion bilinears are ipl(t + e/2, x), ^ 2 (t—e/2, x), which are the same with those 
of ip\ and ipi. Taking into account the time arguments, the following effective Hamiltonian 
Hi +2 



tTNR 
-"kin 



IXLONR I 
-"FV "I" ' 



is derived, where the kinetic term is 



-"kin 



, M \ ~ + / M 

d 3 xZ M l^H Ul + *!>t\ | 02 

-M / -M 



2M 







^~ 



+ y d 3 x^ |^ 

and the interaction Hamiltonian Hpy NR is given by 

Hft™ = H {u) + H {22) + H {12) , 
z 2 r 

H(n) = -f J d 3 xd 3 yV(x - y)T sym {jf(t,x)jf(t,y)} 



^ 2M 



^ 

2M u 



(35) 



// 



(22) 



// 



(12) 



y rfW y v*(f-f)T sym {jf(t,f)jf(t,f)} 

-zZ 9 2 | d 3 xd 3 yD(x - y)T sym {jf(t, x)jf(t, y)} 
Z 2 



(36) 
(37) 

(38) 
(39) 



+ _X J d 3 xd 3 yV x D(x-y)-T syra {jl NR (t,x)jf(t,y)+jf(t,^ 



The time-ordered currents are: 

Tjf = 0^0 1; Tj" = 0} 





fi o y 











of) 

2 



1 
-1 




(40) 
(41) 



Here 0^, "0, M, and g are renormalized fields and parameters and Zm, Z^j, and Z g are 
determined by the renormalization conditions. To be precise, the interaction Hamiltonian 
Eq. (|36|) corresponds to the ladder approximation in the Bethe-Salpeter equation [621.164 165]. 
Cross ladder contributions are of higher order in g 2 and thus can be ignored at leading order 



in the perturbative expansion 68]. 

In the functional Schrodinger equation to be discussed in Sec. IIVI the time-ordered 
effective Hamiltonian H 1+2 provides a convenient way to evolve a state \^(t)) in a 
basis (-01, ^l\^(t)) from an initial condition ^> 2 \^(— oo)) = (0{|ps(— oo)|0|) = 
(0}|ps( — oo)|02)- 111 order to obtain the time evolution of \^(t)) in another basis 
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(Qi(c)iQ2(c)\^(t)) = ( ( 5i(c)IPs , (^)|Q2(c))' is convenient to reorder H 1+2 into a normal- 



ordered product for Q$ c \ and Q^' 



B. Renormalization 

Since the operators in the effective Hamiltonian are time-ordered with respect to infinites- 
imal time differences, interaction terms are expanded in normal order using Wick's theorem 
with the following non-zero contractions: 



: ^ia(x)^(y) ■■ = {■ if>2a(x)ipip(y) : j x (-1) 

0(x - Vo) 



S(x-y). (42) 

ad 



-0(y - x 

The minus sign for : ijj2a(x)'4>2i3(y) ■ d er i ves from the canonical commutation relation: 

{if) 2 a(t,x),i>lp(t,y)} = {$l a (t,x),i> 2 p(t,y)} = -8 a p8{x-y). (43) 

Due to the matrix structure of the Wick contraction, one only obtains a finite contribution if 
the terms are either time-ordered or anti time-ordered the contractions must be time-ordered 
or anti-time-ordered in order to get finite contribution. Therefore any contraction with a 
fermion loop vanishes, which results from taking the quenched approximation. Any contrac- 
tion of a current operator also vanishes because of tr t a = 0. (In QED, this contribution is 
infinitely large because of the filled Dirac sea but is canceled by the background charge.) 
The normal-ordered product of the kinetic term is 

H Z = J (PxZuM {qJQ, + Q\ c Q lc + QlQ 2 + QL&c} 

+ i/:*)* +0 M~£)** (44) 
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and those for the interaction terms are 
Z 2 



H {ll) = -fj d 3 xd 3 yV(x-y)N{jf(x)jf(y)} 
Z 2 C f 

+ / d 3 xV(0)(Q\Q 1 + Q\ c Q lc ) + const., (45) 

H {22) = -Sl J d 3 xd 3 yV*(x-y)N{jf(x)jf(y)} 

z 2 c r 

+ / d 3 xV*(0)(QlQ 2 + Q\ c Q 2c ) + const., (46) 

#(i2) = -iZ] J d 3 xd 3 yD(x-y)N{jf(x)jf(y)} 

+ ^ J d 3 xd 3 yV x D(x-y)-N{jl NR (x)jf(y)+jf^ (47) 



where N denotes the operation of taking a normal-ordered product and we use C-p = (N% — 
1)/2N C = 4/3. The complex potential at the origin ^(0) consists of a vacuum part and a 
finite-temperature part ^(0) = V vac + V mc d- V vac is the vacuum electric potential energy at 
short distance, which is real and divergent, while V me d is the finite-temperature contribution 
to the potential, which is in general complex and finite. 

Since the typical momentum scale of heavy quarks in thermal medium is around a/ MT 
and the typical momentum exchange with medium particles happens at around T, we set 



the renormalization scale at k ~ yMT for the heavy quark kinetic term and at k ~ T for 
the interaction term. First let us renormalize the potential. Actually there is no coupling 
renormalization since the heavy quark internal lines stem from terms we have neglected: 
the many-body interactions with more than 2 heavy quarks involved and the cross ladder 
contributions consisting of 2-body interaction. Therefore we need to choose Z g (T) so that 
it normalizes the vacuum potential V(r) at T = properly at the renormalization scale 
r ~ 1/T. At leading order in the perturbative expansion, we can take Z g (T) = gi^/g^g 
or equivalently gbare = g{T) with Z g (T) = 1. The renormalization condition for the kinetic 



term is that at k ~ a/ MT the effective Hamiltonian describes the non-relativistic heavy 
quark with current mass M at that scale, which is possible only when M 3> T. To be 
explicit, 

Z,(T) = 1, Z M (T) + UT) l^ V - = 1- (48) 
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// 



1+2 



drx 



+ / d 3 x 



The renormalized effective Hamiltonian is finally obtained as 

\m (QlQi + Ql c Qi c ) 

a*M (q\Q 2 + Q 2c Q 2c ) 



(49) 



+ 2 
+ ■ 



d 3 xd 3 y 



v(x - y)N {jf(x)jf(y}} -v*(x- y)N {jf(x)jf (y)} 



-2^(x-y)N{jf(f)jf(y)} 
-±V x d(x-$) -n{^ 



N R (^)j 2 a °( 



;o0| 



with a = 1 + Z 2 C F V med /2M (a = 1 + Z e 2 ^ med /2M for QED), v(r) = Z 2 V(r) } and d(r) = 
Z 2 D{f). The calculation of the gluonic two-point functions G^ (lj, k) and Gq (w, /c) at 
leading order in the hard thermal loop resummed perturbation theory 66( yields (adopting 
coupling renormalization (/bare = g(T) and Z g {T) = 1) 

g(T) 2 [e-^ r „ Z" 00 „ sinfcr 



47T 

d(r) = Im[w(r)] 



- - 2iTuu 
9{T?ulT 



dk 



r(k 2 + ul) 2 



C F g{Tf 



2tt 



-oj d - iT), 



dk 



sin kr 



r(k 2 



2\2> 



2M 4vr 

where ul = g(T) 2 T 2 (N c + N { /2)/3 for QCD and ul = e(T) 2 T 2 /3 for QED. 



(50) 
(51) 
(52) 



IV. REAL-TIME DYNAMICS 



In this section, we derive the master equation for the reduced density matrix of heavy 
quarks and also the time-evolution equation for the forward correlator from the renormalized 
effective Hamiltonian Hi +2 obtained in the last section. A stochastic description equivalent 
to the real-time dynamics is also provided. 
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A. Functional Schrodinger equation for the reduced density matrix 



First let us start with a reduced density matrix p$(t) in the Schrodinger picture: 



p s (t) = Tr E {C/(t,-cx))p tot t/ t (t,-oo)}, 
U(t 1 ,t 2 ) = exp {-i(tx - t 2 )HQCT>} , 



(53) 
(54) 



where p tot = p^ q ® Ps is the initial density matrix for the total system and U{ti,t 2 ) is a 
time-evolution operator with the QCD Hamiltonian Hqcb- Note that the reduced density 
matrix operates only in the heavy quark Fock space. The matrix element of ps(t) between 
the heavy quark coherent states (Ql, Q\ c \ and \Q 2 , Q^) is given as a functional of the heavy 
quark fields by 



Ps 



t, Qi, Qic Q 2 , Q 2c 
(Ql, Qu 



(Q*i,Ql c \Ps(t)\Q* 2 ,Q* 2c ) = (QhQl c ,Q* 2 ,Q 2c Mt)), 



(Q\ exp 



- / d 3 x{Q(x)Ql(x) + Q c (x)Ql(x)} 



(55) 
(56) 



IQIiQzc) = ex P 



■ / d 3 x\Q;(x)Q^x) + QUx)Ql(x] 



\n). (57) 



Here \ fl) is the heavy quark vacuum state. The time-evolution of ps t, Ql, Q\ c , Q2, Q 2c 
obtained by the functional Schrodinger equation for (Ql,Qi c ,Q 2 ,Q 2c \^)'- 



is 



d 

% Ft ps 



t, Qi, Qi c , Q 2 , Q 2c 



H\ +2 

xps 



8 8 

Que) ~^ Q*(c),Qi(c) -> An * , Q 2 (c) ~^ Qt(c)7 Q2(c) -» - ~ 
t, Qi, Qi c i Q 2 i Q 2c 



(58) 



B. Master equation for the reduced density matrix 



By functionally differentiating p s t, Q\, Q\ c , Q^, Q\ c , we obtain, for example, the reduced 
density matrix for a single heavy quark system: 



p%(t,x,y) = (x,i\p Q (t)\y,j) oc (n\Q\x)p s (t)Q^(y)\n) 
5 5 

-Ps 



SQT(x)5Qi*(y)' 



t, Qi, Qi c , Q 2 i Q 



2c 



(59) 
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and that for a heavy quark-antiquark system: 
/°gg e (*>^i>32,j/i,&) = (xi:i;x 2 J\pQ(t)\yi,k;y 2 ,l) 

5 5 5 5 

-ps 



(60) 



5Q\*{x l ) 5QH(x 2 ) 5Q l * c (y 2 ) 8®?{&) 



t, Qi, Qi c , Q 2 , Q 2c 



The master equation for each of these reduced density matrices is therefore obtained from 



the functional Schrodinger equation for p$ t,Q\,Q\ c ,Q2,Q 2c in Eq. ( 1581) . 

For example, let us examine the real-time dynamics of a single heavy quark system at 
high temperature. The master equation for a single heavy quark system reads 
. d 



'at' 



(t,x,y) = < (a - a*)M + 



V 2 - V 2 

x v y 



2M 



, 1 fx X 

+- I 5ij5 kl - -j^- 



-id(x — y) + 



%(t,x,y) (61) 
V x d(x-y) V x - V y 1 kl 



AT 



iM 



and that for the color-traced reduced density matrix pQ(t,x,y) = pQ(t,x,y) is given by 



9 , - 
i-Qj.PQ\t,x,y) 



a-a*)M + 
Cp { —id(x — y) + 



V 2 - V 2 

v x v y 



PQ(t,x,y) 



2M 

V x d{x-y) V x -V v 



,(t,x,y). (62) 



AT iM 

It is interesting to note that taking short distance limit d(x) ~ d(0) + 7x 2 /2 in the master 
equation Eq. (I6"2"j) reduces it to the well-known master equation of the Caldeira-Leggett model 
for the quantum Brownian motion 5l|. The unitarity of the master equation Eq. (162]) is 
easily confirmed using Vd(x)\ x=0 = 0: 

Tr x p Q (t) = J d 3 x(x\p Q (t)\x) = J d 3 xp Q (t,x,x), (63) 

i^^xpqit) = J d 3 xd 3 y5(x - y) \ i^p Q {t,x,y)^ 

= {(a- a*)M - iC F d(6)} Tr x p Q (t) = 0. (64) 

Using the master equation Eq. ( I6"2"j) . the following Ehrenfest relations for averaged expecta- 
tion values (0)(t) = Tr x {pQ(t)0} = J d 3 xd 3 y(x\pQ(t)\y)(y\0\x) are derived: 



4® 



M 



7 



:(P), 



dt x 1 2MT 
±/j?_\ _ SL (/p^_\_3T\ 
dt\2M/ ' MT\\2M/ 2 J 



(65) 
(66) 
(67) 
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To be explicit, 

1 



7 = ^V 2 d(x)\ x=0 



C F g(T) 



2 



^ 2( ^aa,Oo(^)U=0 



9 

C F g(T) 2 f d 3 k 



9 



/d 6 h 
j^y 3 k 2 G> afi0 (u = 0,k), (68) 



which is consistent with the leading-order perturbative calculation of the drag force [llHl3|. 
In this way, we can derive the consequences of classical Langevin dynamics through the 
quantum Ehrenfest relations. 

C. Time evolution of the forward correlator 

The forward correlator of a single heavy quark in the Heisenberg picture is defined as 

G>.(t,f) = Tr{e- pH * CD Q i (t,x)Qtt(0,d)} /Tr { e -^ CD } 

« (n\Tr E {p E Q\t,x)Q>\0,d)}\n), (69) 

where the second equation is a good approximation when e~ M l T < 1. In the Schrodinger 
picture, the forward correlator is 

G^(t,x) » ^e-^(n\(n\U\t,0)Q\x)U(t,0)Q^(6)\n)\n)/j2^ En 

n n 

= (n\Q l (x)Tr E {U(t,0)(p E ® Ps )Ui(t,0)}\n), (70) 

p s = Q j \o)\n)(n\. (7i) 

Here we formally take ps = Q^(0)\Q)(Q\ despite the fact that ps cannot be understood as 
a density matrix for any heavy quark system. This formal treatment is applicable since we 
did not use any property of the reduced density matrix in deriving the Hamiltonian Hi +2 . 
By functionally differentiation, we obtain 

GlAt,Z) ~ T7^PskQl,Ql c ,QZ,QZ c ] ^ , (72) 



5Q?(x) 



>l(c)—*2(c)-° 
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with an initial condition Gq 4 (i = 0, x) — 5ij5(x). Similarly, we obtain the forward correlator 
for the quarkonium state 



G> Qcjij (t, x , y; = Ti{e^ H ^Q l (t,x)Qi(t,y)Q l K0,-mQ k K0,r/2)} /Tr{e-^°} 

(73) 



(n\TT E {p E Q\t,x)Qi(t,^Q l ^0,-f/2)Q k \0,f/2)}\Q) 
6 6 

i, Qi, Qi c i Q21 Q2C 



SQf(x) 8Qil( 



-ps 



with an initial condition Gq Qc ■ ■(£ = 0, x, y; r) = 5ikSji5(x — f/2)S(y + r/2). 

The time-evolution equation for the forward correlator is obtained from that of 
Ps t,Ql,Ql c ,Ql,Ql c in Eq. ( 158]) . For example, the time evolution of forward correlators 



of single heavy quark and quarkonium states is given by 



d 



i— GQ Qc ,ij(t,x,y;r) 



d_ 

dt 



2M 

V 2 . + V 2 



2M 
2 (M« 



aM)G> ti (t,x) 



(74) 



' +2aM)G% Qctij {t,x,y;r) 



IvT 



v{x-y)G^ Q Jt,x,y;r). 



(75) 



Projecting the heavy quark state on the color singlet GQQ cl {t,x,y;r) = GQQ cii {t,x,y;r), 
we get 



9 



V 2 

V ... 



2M 



+ 2aM - C F v{x-y) } G>(t,x,y;r), (76) 



from which the complex color-singlet potential can be read off as t>i(r) = 2{a — l)M — C F v{r). 
Using the perturbative results in Eqs. fl50|) . ( 15~T|) . and (152]) . the complex color-singlet potential 
at leading order in the perturbative expansion is 



vi(r) 



C F g(T) 
Air 



2 r 



-o) D r 



w D + 



zT <! 1 - 2u;£ 



00 sinfcr 
dk 



o ">(F + u; 2 ) 2 



(77) 



which agrees with the complex potential obtained by the spectral decomposition of the real- 



time thermal Wilson loop in leading-order perturbation theory 28|. Projecting onto the 

(6ij/N c )C 

v(x - y) 



color-octet state G% QcjS (t, x, y; f) = G^ Qc ij (t, x, y; f) - (5ij /N c )G> Qc kk (t, x, y; f), we get 



d 



v^ + v 2 

" V2aM + 



G^ 8 (t,x,y;r). (78) 



2M 2N C 
The complex color-octet potential Vs(r) = 2(a — 1)M + v(r)/2N c is therefore obtained as 



v 8 {r) 



Air 



-C F u D + — — %T 

2i\ r r 



D sin At 
dk 



'r(k 2 + ul) 2 



(79) 
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in leading-order perturbation theory. It is interesting to note that the imaginary parts of 



complex color-singlet and color-octet potentials are different at the origin (Im t>i(0) 







and Im 



v s (0) 



—T(N c /2)g(T) 2 /An) while they are the same at infinitely long distance 
Im[f!(oo)] = Im[u 8 (oo)] = — TC-pg{T) 2 /Ait. Since the imaginary part reflects Landau 
damping due to the collisions between the heavy quark and the medium gluons and light 
quarks, it is natural to obtain a vanishing (finite) imaginary part in the complex color-singlet 
(color-octet) potential at the origin. 



D. Stochastic description 



The effective action S1+2 can be expressed partly as a stochastic average using both real 
(£"(x)) and complex (£" (x), £f (x)) Gaussian white noise [691 ]: 



+2 



e 

,• cLONR 



g kin + l °FV + 



exp 



- / dt I d 3 xd 3 yRe[V(x-y)} (jf(t, x)jf (t, y) - jf(t, x)jf (t, y)) 



(80) 
(81) 



x ( exp 



i I d x 



mx)+tf(x)}jf(x) - £$(x)V • jU r (x 



where all the finite noise correlations are 



-5 ab 5(t- s)D(x-y), 
-<5<%5(t- s)2D{x -y). 



(82) 
(83) 



Note that D{f) is negative definite so that the probability density for the noise can be 
normalized. 

Following the same procedure as in Sec. IIIH the normal-ordered renormalized effective 
stochastic Hamiltonian H® 2 * s obtained as a sum of H[® and H!p which consist only of 
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Qim and Q2L) respectively: 



e kin +U FV 



HI 



(0 



^exp 



1 At H i 



K) 
+2 



1 



+^ I d 3 xd 3 y 



zr(0 
-"1+2 



Re [a] M (qIQx + Q f lc Q lc ) + [q\ (-g) Qi + Q f lc (- § ) Q 1( 
+ {£(*, 2) + «)} N {jf (£)} - ^f(t, 2) V • N (j? )NRl 

ReKx-$]N{jf(£)if($} 



34) 



•r, 



35) 



if 



(0 



Re[a]M (q|Q 2 + QLQ 



2c 



Q 



vi 

2M 



Q2 + Q 1 " 



2c 



2Af 



Q 



2c 



- {£(*, x) + £?(t, £)} N {jf (£)} - j±,£r (*, £) V • N jl 



2,NR 



d 3 xd 3 y 



Re[v(x-y)]N{jf(x)jf(y)} 
with renormalized noise correlations 

(C(t,x)e v (s,y)) = -5 ab 5 ts d(x-y)/At, 



36) 



The time evolution of p s Qi, Qi c ; Q25 Q 



2c 



is calculated by 



Ps 



t + At,Ql,Q* lc ,Q* 2 ,Q 



2c 



iAtH 



(€> 

1 + 2 



Ps 



t, Ql, Qlc; ^2c 



(87) 
(88) 

(89) 



iAtH\ 



(0 
+2 



At 2 



(0 



Ps 



it, Qi , Qi C ! Q2' ^2c 



Note that stochastic variables are of order 0(1/ At 1 / 2 ) and therefore terms up to iij+2 are 



necessary even in the limit At — > 0. Since H\ 
time-evolution operator is given by 



(?) 
+2 



(0 



ifP and 



' ' ? -"2 



0, the 



iAtff} 



(5) 
+2 



At 5 



/A/// 



(0 



At 2 



(0 



At 2 



// 



in the limit At — ?■ [70j . Here only the real noise correlation remains finite in ( H{ 



r(0' 



(90) 
and 



// 



. This decomposition enables us to describe the time evolution of the system in 



terms of a stochastic process with real and complex noises. 

For example, the stochastic evolution of a single heavy quark system is obtained as 

r(0 (z\TTtt) f.t\\ M, 
G%, k (t,x) 



p^t + At,x,y) = {U^U^^)^^ 
G> ti (t + At,x) = (u { Q %(x) 



(91) 
(92) 
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with U Q&(x) and U S!ik(x) defined by 



(0 



6Q?(x) 
5 

WW) 



1 - iAtH\ 



(0 



1 - zAtK 



(0 



At 2 
IT 

At 2 



(O 1 



H. 



J 8Q l *{x) 



(93) 
(94) 



-2(c) 



Here and H^p are functional operators obtained from H[ and ff^ by the same 



substitution as in Eq. ( 158]) . Explicit matrix forms are 

V 2 



r^(2) = i-iAt 



^(f) = 1 + zAt 



aM 



2M 
" 2M 



te a (x)+^(x)}f+ V^(x) f 



8MT 



{0) + C(x)}t a *+ Vffi * 



8MT 



(95) 
(96) 



Using the time-evolution operators Ug\ k (x) and Uq^x), the following stochastic 
Schrodinger equations are obtained: 



¥ Q (t + At,x) = U Q %(x)* k Q (t,x), 
%(t + At,x) = U { Q %(x)¥ Q (t,x), 



(97) 
(98) 



with which the density matrix PQ(t,x,y) and the forward correlator Gn^t, x) are given by 

(99) 
(100) 



p%(t,x,y) = ^ Q (t,x)%(t,y) 
G>^t,x) oc (%(t,x)) r 



Using this stochastic description and taking the classical limit, we can derive the Langevin 
equation, which we summarize in the Appendix [B] A stochastic description for systems with 
arbitrary number of heavy quarks can also be obtained in a similar manner: for a single 
heavy antiquark system, Uq^(x) {Uq^{x)) can be obtained by substituting — t a * for t a (t a 
for —t a *) in Uq\x) (Uq\x)); for a heavy quark-antiquark system UqQ c (x,y) = Uq\x) + 
U®{#) - 1 - iAtv{x-y){t a )Q ® H°*)q c and U^JS, y) = U®(x) + U { £(y) - 1 + iAtv*(x - 
y)(-t a *)Q®(t a ) Qc - 

It is interesting to observe here that Uq^x) ^ Uq^*Ax) and &Q(t,x) ^ \l/g(t,x). When 
M = oo, the complex noises £f(x) and £f (x) can be eliminated so that Uq^-{x) = Uq^*(x) 
and ^ l q(t,x) = ^q{t, x) hold. This is a purely stochastic process with Uf{x) = 
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exp[— iAt£f(x)t a ], which is not an irreversible process. Since the complex noises and the 
irreversible drag process originate from the same term, i.e. ^7inrJ2° + Ji°J2NrJ m ^ ne 
fective action 5 , i +2 , it is natural to obtain independent stochastic time-evolution equations 
for ty l Q(t,x) and x) in the presence of the complex noises. This gives us a theoret- 

ical foundation for the stochastic description of Ref. 50| and also illustrates how it must 
be extended to the finite heavy quark mass case. In summary, stochastic processes with 
the real noise encode only the fluctuations while those with complex noises incorporate the 
irreversible processes caused by the drag force in the classical Langevin dynamics. 

V. CONCLUSION AND OUTLOOK 

In this paper, we have studied the real-time dynamics of heavy quarks as an open quantum 
system in thermal QCD medium. We base our approach on the closed-time path formalism of 
non-equilibrium quantum field theory and have worked in leading-order perturbation theory 
and in the non-relativistic limit to make calculation tractable. We have shown that all 
time-evolution equations can be derived from the renormalized effective Hamiltonian -H1+2 
for doubled degrees of freedom of heavy quark fields: the master equation for the reduced 
density matrix Pqit, x, y), the time-evolution equation of forward correlator Gq Qc ^(t, x, y; f), 
and the stochastic Schrodinger equation for forward-propagating {^fqit, x)) and backward- 
propagating (^/q{t,x)) wave functions. The renormalized effective Hamiltonian governs the 
time-evolution of the functional reduced density matrix p$ t, Ql, Ql c , Q%, Q2C i n the form 
of a functional Schrodinger equation. In leading-order perturbation theory, the renormalized 
effective Hamiltonian depends only on two real functions, i.e. the real and imaginary parts 
of the complex potential v(f). Note that the imaginary part of the complex potential also 
determines the drag force in perturbation theory. 

In the stochastic description, the real part of the complex potential Re [v{r)\ represents 
the interaction between heavy quarks while the imaginary part of the complex potential 
d(r) = lm[v(r)} determines the real noise correlation. It is also worth noting that the 
stochastic description introduces complex noises, whose coupling to heavy quark fields is 
not uniquely determined. This indicates that the complex noises are not actual physical 
observables but only represent mathematically the friction process in the stochastic decom- 
position. From a practical point of view, however, introducing complex noises enables us 
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to numerically simulate heavy quark systems, for example a single heavy quark system, by 
stochastically evolving two 3-dimensional functions ^^(t, af) and ^^(t, af). Compared to 
numerical simulation of a 6-dimensional density matrix Pq(t, x, y), the use of complex noises 
might thus be advantageous. 

Although our derivation relies on perturbation theory, it should be helpful to construct 
phenomenological and effective description in the non-perturbative region, perhaps by mod- 
eling the renormalized effective Hamiltonian -H1+2, from which the various time-evolution 
equations are derived. One strategy will be to expand JT1+2 in terms of heavy quark bilin- 
ear Q\( c sQi( c ) an d Q\^Q2(c) with each coefficient to be determined by matching. Another 
strategy may be to use the reduced density matrix in thermal equilibrium, which is a static 
solution of the master equation. Some information must be gained by comparing the thermo- 
dynamic quantities, e.g. the Polyakov loop and its correlator, calculated non-perturbatively 
and those obtained from the reduced density matrix in equilibrium. Our analyses have 
shown that thermal fluctuations and the friction process are both indispensable ingredients 
to the real-time dynamics of heavy quarks and that the number conservations of heavy quark 
and heavy antiquark do strictly hold in spite of the complex potential. In constructing a 
model, it should also possess these properties. 

In order to fully appreciate the non-perturbative dynamics of the deconfinement transition 
through heavy quark systems, there still exists a following nontrivial challenge. How should 
the dynamics of heavy quarks be described in the confined phase, where heavy quarks are 
bound not only in quarkonium states but also in heavy-light mesons? More specifically, how 
is our potential-based description connected to the state-based description by the quantum 
optical master equation. This seems to be very challenging especially when a heavy quark is 
bound in a heavy-light meson because there is no potential between heavy quarks and light 
quarks in our description. It might be necessary to include light quarks that are resonantly 
interacting with heavy quarks into the open quantum systems. 
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Appendix A: Relations among the two-point functions 

Here we derive some of the relations among the gluonic two-point functions in the in- 
stantaneous approximation. From the definition, it follows that 

<5afe,oo( f -y) = -iG^ b>00 (x-y) + Gf b ^(x-y), (Al) 
<5a6,oo (x-V) = i&abfio (z - y) + G> b4iu (x - y) , ( A2) 

Gab,Oo(x-V) = G ab*nv(x-y)- (A3) 

Using the Kubo-Martin-Schwinger (KMS) condition G^ b00 (t, x — y) = G^ b00 (t — if3,x — y), 
analyticity of the forward correlator G> b00 (t, x — y) in a stripe —ij3 < x° — y° < in the 
complex time plane and assuming that the damping of the correlation at \x° — y°\ — > oo, 
the instantaneous approximation of the retarded propagator yields an Euclidean correlator: 

/oo 
dx°9(x° - y°){[A%(x), A b (y)]} 
-oo 

dr{A a (t = -ir,x)A b (t = 0,y)). (A4) 

Obviously G^ b00 (x — y) is real. Using the KMS condition, we can derive G> b00 (x — y) = 
G< b00 (x — y) and thus G^ boo (x — y) is also real. The above results give the relations among 
Gab i00 (x- y)s. 

The spectral function of gluon fields, which is a real and odd function of u, is expressed 
in terms of the forward correlator and retarded propagator: 

/oo 
dx°e^ a -^([A a Q {x),Al{y)}) 
-oo 

= (l-e-^)G> fi0 (u;,x-y) 
= (e^-l)G< bj00 (u,x-y) 

= 2lm[G% (u,x-y)]. (A5) 
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Using the fact that G'^ 0Q (x — y) is purely imaginary, which derives from the definition of 
G ab,oo( u ' x-y),we obtain 







G ab,oo(%-y) = T—a ab)00 (uj, x — y) 



oj=0 
/R 



G% m {x - y) = -G"i 00 (f -y) = -iG' a u 00 (x - y) 



Id _ 
2duj (Xab ' 00 ^ UJ,X ~ V ' 



LU = 



(A6) 



(A7) 



which completes the derivation of all the relations among two-point functions. 



Appendix B: Derivation of the classical Langevin equation 

We show here how the classical Langevin equation is derived from the stochastic descrip- 
tion. For simplicity, we derive the Langevin equation for muons in a QED plasma composed 
of photons and electrons. Let us start from the stochastic representation of the master 
equation Eq. ( )9T|) adapted to QED: 



p Q (t + At,x,y) = (u%\x)U®(tf)^p Q (t,x,&). 



(Bl) 



Using the stochastic evolution operators: 

V 2 



U { i\x) 



U^(x) 



1 -iAt 



1 + iAt 



aM 



V 



2M ' * v " I " v ' ) 8MT 
V 2 „ „ r^, ,! V 



a * M ~ 2M + + + 



8MT 



(B2) 
(B3) 



and taking the noise average only for the complex noises and &(x), the master equation 
is derived as 



p Q {t + At, x, y) - p Q (t, x, y) 

v?-v* 

-iAt , 

-id(x -y) + 



+ Zt(x) - Uy) 



(a - a*)M 2M 

V x d{x-y) Vk-V 



pQ(t,x,y). 



(B4) 



iAI 



Regarding this stochastic master equation: 



, (t + At , x , y\ £ r ) - p Q (t , x , y; £ r 

v?-v 2 

la — a" )M — 

-iAt 



2M - + Ux)-Uy) 

id(x-y) + ^^-^ 



(B5) 



26 



with pQ(t,x,y) = (pQ(t,x,y;^ r ))^, the following stochastic Ehrenfest relations are derived 
for each noise history: 

(S)(f + Af;&) = (^(t^J + At^M), (B6) 
(p)(t + At;& = (p)(t;^)-At^;(p)(t;6)-At J d 3 xV^-p Q {t,x,x;^). (B7) 

In the classical limit, or more precisely, when the probability density of the muon wave 
packet is localized around x ~ (x)(t; £ r ) with finite spatial size over which the noise £, r {x) is 
effectively constant, we can make an approximation x, x; £ r ) ~ 5(x — (x)(t; £ r )). Finally 
we obtain the classical Langevin equation for momentum diffusion, 

<#(f + Af;&) = (p)(t^ v )-At^(m^r) + Atf(x = (x)(t-,^)), (B8) 
/(f) = -Vf r (f), (B9) 
{fi(x)fj(*)) = --fa i Q-d(x-$)\ g= t/At = 1 5 ij /At. (BIO) 
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since VD(x) ~ T ■ D(x) while VQ^(x) ~ ^MT ■ Q?l(x). 

[68] Such contributions can be included systematically by the correspondence: exp [— ii3" 1+2 t] = 
T(sym) exp —i Jq dt' (i^n + -tfpy NR + • • • ) , where each order of the expansion of the right 
hand side can produce terms proportional to t. Thus even if we take into account terms only 
up to -ffpy NR , this expansion gives iJi+2 as an expansion in terms of g 2 , whose leading order 
(higher order) corresponds to the ladder (cross ladder) contribution to -Hi+2- 

[69] This stochastic expression is not unique. There is a freedom in how the complex noises couple 
to the heavy quark fields. For instance, we can also take c • (,fjf° and c" 1 • £"*V • j% NR instead 
of £fji° and £"*V • J2NR w ith some complex coefficient c 7^ 0. 

[70] In the stochastic description, the ladder approximation corresponds to taking At — > in the 
expansion of e- iAtH &. The cross ladder contributions, which turn out to be of higher order 
in g, are recovered by taking into account the finite correlation time At 7^ in the expansion, 
which is not necessarily an expansion in terms of small dimensionless number. See 
similar argument. 
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